% 创建时间：2025-06-20 22:38:43
% 作者：兪馳
clc;
clear;

A=[1 , 0 , 4 , 0; 
    0 , -3, 0 ,0; 
    -2 , 2 , -3,0; 
    0 ,0 ,0, -2]; 
[V,D]=eig(sym(A))

disp('----------------------------------------')

E=eye(size(V));
k=4;
E(:,k)=E(:,k(end:-1:1));%E(:,k)表示选择矩阵E的所有行和第k列；E(:,k(end:-1:1)) 则对列索引进行反转操作，即从k的最后一个元素开始，到第一个元素结束。
V=V*E;
J=inv(V)*A*V%这里在进行对角化
A=[1,0,4,0;
    0,-3,0,0;
    -2,2,-3,0;
    0,0,0,-2];
[V,D]=eig(sym(A));
[V1,A1]=jordan_real(sym(A))